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ABSTRACT 

The thermal evolution of neutron stars is coupled to their spin down and the resulting 
changes in structure and chemical composition. This coupling correlates stellar surface 
temperatures with rotational state as well as time. We report an extensive investigation 
of the coupling between spin down and cooling for hybrid stars which undergo a phase 
transition to deconfined quark matter at the high densities present in stars at low ro- 
tation frequencies. The thermal balance of neutron stars is re-analyzed to incorporate 
phase transitions and the related latent heat self-consistently, and numerical calcula- 
tions are undertaken to simultaneously evolve the stellar structure and temperature 
distribution. We find that the changes in stellar structure and chemical composition 
with the introduction of a pure quark matter phase in the core delay the cooling and 
produce a period of increasing surface temperature for strongly superfluid stars of strong 
and intermediate magnetic field strength. The latent heat of deconfinement is found to 
reinforce this signature if quark matter is superfluid and it can dominate the thermal 
balance during the formation of a pure quark matter core. At other times it is less 
important and does not significantly change the thermal evolution. 

Subject headings: Stars:neutron — stars:rotation — dense matter — equation of state 
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1. Introduction 

The chemical composition of neutron stars 
at densities beyond nuclear saturation re- 
mains uncertain with alternatives ranging 
from purely nucleonic compositions through 
hyperon or meson conden sates to deconfined 
quark matter - see e.g. IWeber and 
Page Reddvl (|2006l l for recent reviews with 
emphasis on quark matter. A future un- 
derstanding of neutron star structure gained 
through confrontation of theoretical mod- 
els with the now steadily growing body of 
observational facts will therefore simultane- 
ously constrain fundame ntal elements of par- 
ticle and nuclear physics (jLattimer &: Prakash 
20071 ). The interlinked processes of spin 
down and thermal cooling present intriguing 
prospects of gaining insight in the properties 
of matter in neutron star cores by confronta- 
tion with soft X-ray observations of thermal 
radiation from neutron star surfaces as they 
both depend sensitively on and to some extent 
determine the chemical composition. 

As neutron stars spin down and contract, 
their structure and composition change with 
the increasing density - drastically if phase 
boundaries are crossed so new forms of mat- 
ter become possible. We shall here investigate 
how this influences the thermal evolution of 
hybrid stars which contain large amounts of 
deconfined quark matter. The increasing den- 
sity and changing chemical composition fur- 
ther imply additional entropy production in 
bulk and the release of latent heat as particles 
cross any phase boundaries present. We there- 
fore re-analyze the thermal equilibrium of 
compact stars to show how mixed phases may 
be incorporated. We thus arrive at a natural 
description of the latent heat of phase transi- 
tions in compact stars, but also find - through 
direct numerical calculations - that unless the 
stellar structure changes very rapidly the ef- 
fects of latent heat on the thermal evolution 



are insignificant when compared to those of 
the changing chemical composition and the 
surface area reduction. 

Neutron stars are extremely compact ob- 
jects and densities in their cores reach well in 
excess of the nuclear saturation density. At 
such densities the distance between particles 
is on the order of the characteristic range of 
the nuclear forces. Therefore, as was stressed 
recently by Baym ( 200?! ). perturbative treat- 
ments in terms of few- or even many-body 
forces - although highly successful in describ- 
ing the properties of matter below nuclear sat- 
uration - are no longer well defined in the core 
of neutron stars. Further, the relevant degrees 
of freedom should include the appearance of 
hyperons and possibly deconfined quarks, so 
treatments in terms of nucleons alone must 
also be seen as approximative. Hyperons are 
expected to appear at densities around 2—3po, 
where po = 0.153 fm~^ is the baryon density 
at nuclear saturation. The appearance of hy- 
perons so softens the equation of state that 
purely hadronic equations of state may not 
allow stable models compatible with the accu- 
rately measur ed masses of neutron stars in bi- 
nary sv stems (IBaldo et al. [good : I ISchulze et al 



200i). ISchulze et al.l (|200fil ) further demon- 



strate that this conclusion is highly robust 
with respect to different assumptions about 
hyperon interactions and the nucleonic equa- 
tion of state. At best these studies are strong 
arguments against purely hadronic composi- 
tions and they certainly do show the relevance 
of considering alternatives such transition 
at high densities to deconfined quark matter. 

Quark matter represents an entirely new 
type of matter - as opposed to just an ad- 
ditional degree of freedom as in the case of 
hyperons in the hadronic phase - and it can- 
not be assumed to softe n the equation of 
state to the same extent ( Alford et al. 20081 : 
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Schaffner-BielichI l2007l ) . Hybrid stars can be 
be consistent with the high masses and radii 
indica t ed by recent observ ati ons (e. g. 



Qzell 



\2OOd ). iFreire et al.1 (|2008bl lah: iFreird (|2008l l1 

and the quark matter equation of state can 
fulfin the constraints imposed by heavy-ion 
cohision transverse flow data and pro- 
duction (which nu cleonic equa t ions o f state 
do not). Further, lOrago et aD ( 20081 ) found 
that only stars with a quark matter compo- 
nent can rotate stably without losing angular 
momentum by emission of gravitational waves 
through the r-mode instability at the 1122 Hz 
indicated by recent observatio ns for the X- 
ray t ransient XTE J 1739-285 (|Kaaret et al. 
20071 ). This conclusion is tentative as the 



observed neutron star spin frequency awaits 
confirmation, but it again shows how the com- 
position of neutron stars is an open question 
to be determined by a confrontation of theory 
and observation, and that a deconfined quark 
matter phase remains a viable alternative. 

For definiteness we shall work with the 



equat ion of state suggested by iGlendenning 



(|l992l ): tiereafter the GfgQ equation of state. 



While this equation of state is not sophisti- 
cated in its treatment of the quark matter 
phase, it is illustrative in that it allows a very 
large pure quark matter core with a transi- 
tion through a mixed phase at relatively low 
densities around 2-5po- For our purposes it is 
interesting in that it implies drastic changes in 
composition with spin down - the pure quark 
matter core disappears at high spin frequen- 
cies for instance - and it therefore represents 
something close to a limiting case. We shall 
then be able to investigate both the effects of 
steady conversion of hadronic matter to quark 
matter and the sudden appearance of a pure 
quark matter core in the hadronic phase. 

In the so-called minimal scenario, which 
exclude s exotic and very rapid neutrino pro- 
cesses (jPage et al.l |2004| ). the internal tem- 



perature of neutron stars drops from beyond 
10^*^ K to around 10^ K within a few minutes 
after their birth, and neutrino cooling con- 
tinues to dominate for at least the next few 
thousand years until the internal temperature 
has dropped below 10® K and photon coo l- 
ing takes over (see, e.g. Page et al. ( 20041 ) : 
Yakovlev k Pethickl (|2004l ) for reviews). If 



the highly efficient direct Urea neutrino emis- 
sion (i.e., essentially beta decays, n ^ p + 
e~ + i>e and related processes in quark matter) 
is active, the stars may cool very rapidly and 
reach very low temperatures on a timescale of 
a few hundred years. As we shall see the nu- 
clear direct Urea process is active in the mixed 
phase of hybrid stars, and may thus control 
the thermal evolution. These processes may 
be suppressed by pairing of the participating 
particles however, and further the extent of 
the mixed phase depends strongly on the ro- 
tation frequency thus giving rise to a diverse 
range of possible cooling paths. 

The latent heat of the phase transition 
and the related release of entropy in bulk 
with changing density are generally found 
to be of smaller importance than the chang- 
ing structure and chemical composition. It 
does not significantly delay or enhance the 
cooling, although it does briefly balance or 
even dominate the cooling terms when the 
pure quark matter core is first formed for cer- 
tain choices of stellar parameters. The latent 
heat of the quark matter p hase transition was 



previously considered b y iMiao &: Xiao-Pine 



( 2007 '): lMiao et al.l (|2007l ') and. lXiaoping et al. 
(j2008l ). These authors took a different ap- 



proach to calculate the latent heat than what 
is discussed below and found very significant 
heating terms which we do not recover. 

The work of iReisenegger (I1995I ) and lately 
Fernandez &: Reisenegger I2OO5I I considered 



the related process of roto-chemical heating 
for neutron stars in which weak reactions are 
driven out of equilibrium by the changing den- 
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sity. The resulting release of energy was found 
able to maintain old stars at relatively high 
temperatures determined by their rate of spin 
down. This term naturally appears in our 
equations for the thermal balance and should 
be included in a full model. It is to some ex- 
tent complementary to the effects discussed 
here, but the treatment of weak reactions be- 
yond equilibrium is beyond our scope in this 
work, and we shall assume chemical equilib- 
rium throughout. We return briefly to discuss 
this issue in Sect. [5l 

The link between the thermal evolution of 
neutron stars and their spin down has be- 
come a subject of some interest as an ob- 
servational correlation between the tempera- 
tures and inferred magnetic fields of neutron 



stars was recently discovered by iPons et al 



(|2007l l. This was interpreted as evidence for 
magnetic field decay and detailed work by 
Aguilera et all (|2008l E strongly supports this 
conclusion. An alte rnative interpretation was 
offered recently bv iNiebergal et ah ( 2007 ) in 
terms of magnetic flux expulsion from color- 
flavor locked quark stars (a hypothetical class 
of stars consisting entirely of quark matter 
which in this case is assumed to be absolutely 
stable). Th e prev iously mentioned work of 
Drago et al. ( 20081 ) also indicates a link be- 
tween the spin down and cooling of hybrid 
stars. These correlations - if they can in- 
deed be shown to exist - complement the 
traditional cooling calculations which relate 
only temperature and age, and they may help 
break the degeneracy seen between such calcu- 
lations with different underlying assumptions 
about the state of matter at very high density. 

In the following, we shall first revisit the 
equations of thermal balance for compact 
stars in Sect. [2] to show the effects of a time- 
dependent density and discuss how the pres- 
ence of phase transitions may be included. 



In Sect. [3l we discuss the Gfgg equation of 
state, the resulting stellar models and some 
simple estimates for the additional terms in 
the thermal balance equations in more detail. 
In Sect. [3] we show the results of including 
these terms in spherical isothermal cooling 
calculations. We conclude with a discussion 
in Sect. [5l 

2. Thermal Equilibrium in the Mixed 
Phase 

The thermal evolution of compact stars is 
determined by the equations of local energy 
conservation and transport in the framework 
of general relativity. These equations balance 
any energy radiated away from the star by 
photons and neutrinos against changes in rest 
mass due to nuclear reactions and changes in 
gravitational or internal energy as the stellar 
structure evolves. In the most widely studied 
scenario the stellar structure is assumed con- 
stant and the only energy source available to 
neutron stars is then their original end o wmen t 



of thermal e i iergy (e.g.. I Van Rip eij ( 1991 ): 



20061 1: 



ISchaab et al.' ('1996''): 'Page et al.l (|2004l . 

Yakovlev & Pethick (2004) and references 
therein). Neutrino production in the core and 
photon emission from the stellar surface then 
ensures a monotonical and sometimes very 
rapid cooling of the star depending strongly 
on any superfluid properties of the core. But 
in addition to this a number of powerful en- 
ergy sources may play a role in delaying or 
even reversing the cooling at various stages 



see e.g. 



Schaab et al 



i n a ri eutron stars life _ 

(|l999l ) :br an overview of possible sources and 
their effects. Here we discuss the effects of 
including a phase transition in the energy bal- 
ance, which, as the stellar structure changes, 
must then also include redistribution of en- 
tropy between the two phases as their pro- 
portion changes as well as changes in surface 
and Coulomb energy for transitions through 
a mixed phase. 
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Following iThorne 



(119661) or the equivalent 



discussion in IWebeij (|l999l ) we consider en- 
ergy conservation for a spherical shell inside 
of which is a baryons and which itself con- 
tains 6a baryons. The treatment given in this 
section therefore assumes spherical symmetry 
which is sufficient to show the effects of a 
phase transition on local energy conservation 

- we shall return later to discuss possible con- 
sequences of a multidimensional cooling cal- 
culation. The shell under consideration will 
change its internal energy during a coordinate 
time interval dt by 

o?(internal energy) = 

(amount of rest mass-energy converted to 
internal energy by reactions) 
+ (work done on shell by gravitational 
forces to change its volume during 
quasi-static contraction) 

— (energy radiated, conducted 

or convected away from the shell). (1) 

It is important to note at this point that if 
the two phases are in equilibrium - as should 
certainly be expected for transitions governed 
by strong reactions such as that from hadronic 
to quark matter - there is no binding en- 
ergy involved in the transition and therefore 
no contribution to the first term on the right 
hand side of Eq. ([T]). If this was the case 
and the phase transition did involve a bind- 
ing energy the two phases could not be in 
equilibrium and the star would adjust itself 
accordingly - eventually becoming a strange 
star in the case of the quark matter transition 
if quark matter was assumed bound relative to 
hadronic matter at zero external pressure. In 
hybrid stars this is not the assumption how- 
ever, and so any latent heat evolved or ab- 
sorbed in the phase transition follows from 
the different thermodynamical properties of 
the two phases as with any other phase tran- 
sition. 



In the mixed phase of the Gfgf] equa- 
tion of state regions containing negatively 
charged quark matter appear at densities of 
about 2 times nuclear saturation and dom- 
inate completely at 5 times nuclear satura- 
tion. They allow the hadronic matter to 
lower its isospin asymmetry energy and be- 
come positively charged by including more 
protons with charge neutrality achieved glob- 
ally. The geometry and structure of the 
mixed phase is determined by a balance be- 
tween surface tension and Coulomb repul- 
sion between regions of like charge. For de- 
tails on the phase tr ansition we refer to e.g . 

iHeiselbere et al.' (^199^^ 
Glendenning &: Weber 



Glen denning 
Glend enning 




200n : IVoskresenskv et al.l (|2003l ): lEndo et al. 
20061 ). 

Returning to our spherical shell of baryon 
number 5a we will assume that its volume V 
is large enough to contain a macroscopic num- 
ber of unit cells each containing a region filled 
by the rare phase whose presence may then 
be considered a microscopic property of the 
equation of state. The work done to change 
the shells volume during contraction or expan- 
sion of the star must then include changes in 
surface and Goulomb energy as well as the 
usual pressure term 



dW 



-PdV + adS + dEc 



(2) 



where a is the surface tension, S is the amount 
of surface area in the shell dividing the two 
phases and Eq is the Coulomb energy con- 
tained in the shell. These terms can be sim- 
ilarly included in the first law of thermody- 
namics which we write in the same notation 
(units with G = c = = 1 will be used here 
and throughout this paper) 



de 



p 



dp + Tpds + l^kpdYk 



5a ^ P 5a 



(3) 



where s is the entropy per baryon and 1^ = 
Pk/ pis the fraction of the baryon number den- 
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sity in the form k with k running over ah par- 
ticle species and X^fc ^fc ~ 1- quark 
matter phase each quark contributes by only 
I to the baryon number density and s is then 
the entropy per three quarks. The last two 
terms in Eq. ([3]) are the local densities of sur- 
face and Coulomb energy written in a form 
useful for our purpose. 

Inserting this in Eq. ([T|) and using that the 
amount of energy radiated, conducted or con- 
vected away from the shell can be written in 
terms of the gradient of the total luminosity, 
Ltotj we show in Appendix lAlthat local energy 
balance may be expressed as 



da 



2$\ 



dt 
(4) 

thus giving the contribution to Ltot from the 
shell in terms of the entropy production and 
the difference between the chemical potentials 
of particles participating in any reactions tak- 
ing place in the shell, e^* is the time compo- 
nent of a spherically symmetric metric 

ds^ = -e^^dt" + e^^dr^ + r^dO^ + sin ^^^$2 

(5) 

found from the general relativistic structure 
equations for compact stars. 

Ltot includes the neutrino luminosity, but 
since neutrinos can be taken to immediately 
escape from the star when they are created 
without converting into any other form of en- 
ergy along the way, they fulfill their own sep- 
arate equation of energy conservation 



d 

da " 



(6) 



where is the neutrino luminosity and 
is the neutrino emissivity; the rate per unit 
volume at which neutrino energy is created. 
In neutron stars convection is negligible com- 
pared to electron conduction and photon dif- 
fusion and so the remainder of Ltot can be 



show i i to fulfill a t ransport equation ([Thorne 
19661 : IWebei]ll999l ) 



da 



Le^ 



16cr r3 167r2r4p 



(7) 



where a is the Stefan-Boltzmann constant, k 
is the total thermal conductivity and L = 
Ltot — L^. At the stellar center we must have 
-^^tot(o = 0) = while at the stellar surface L 
must equal the total stellar photon luminos- 
ity which may depend on assumptions about 
properties of the stellar atmosphere or lack 
thereof. 

Eqs. dH)-© with the appropriate bound- 
ary conditions can be solved to evolve the 
thermal structure of a stellar model. They 
have exactly the same form as would be ex- 
pected in the absence of any phase transi- 
tions. However, the phase transition infiu- 
ences the entropy density and this, as we 
shall see, gives rise to additional terms in the 
heat balance including latent heat and the 
surface and Coulomb energies. An equiva- 
lent form of Eq. showing the contribu- 
tions from surface and Coulomb terms more 
explicitly can be found in Eq. (|A6P of Ap- 
pendix [Al Combining Eqs. ^ and ([6]), not- 
ing that Cv = pT{ds/dT)v and assuming con- 
stant structure and chemical composition we 
also recover the standard cooling equation for 
static stars 



-P 



-1 



2$ 



d(re^ 



(8) 



Eq. dH) is useful because it allows a particu- 
larly simple analysis in the presence of a phase 
transition. The first thing to note is, that 
particles crossing the phase boundary would 
not contribute to the second term on the right 
hand side of Eq. @ if the two phases are in 
or close to equilibrium, because their chemi- 
cal potentials (or those of their constituents) 
must then be continuous across the phase 
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boundary. In the following we therefore ne- 
glect this term, but we shall return to discuss 
its potential importance for reactions not in 
equilibrium. The first term is a different mat- 
ter however. The entropy per particle is a 
function of density, temperature and chem- 
ical composition which is not required to be 
continuous across the phase boundary, so par- 
ticles making the transition will release or ab- 
sorb heat accordingly. To see how this works 
we write s as a sum with bulk contributions 
from each phase according to their volume or 
mass fraction as well as contributions from the 
surface and Coulomb energies 

s = -{{l-x.)Si + x.S2 + Ss + Sc) (9) 
P 

= {l-Xp)si+XpS2 + ^{Ss + Sc) (10) 

where p = {1 — Xv)pi + XvP2 is the average of 
the two particle densities pi^2 weighed by the 
volume fraction of the dense phase, Xv Si^2 
is the bulk entropy density in each phase and 
Ss and 5c are the surface and Coulomb con- 
tributions to the entropy density respectively. 
We have further introduced the baryon num- 
ber fraction Xp of baryons in the dense phase 
in a sample of matter related to the volume 
fraction by xp = XviP2/p)i as well as the en- 
tropies per baryon in the respective phases, 
si^2 = 81^2/ Pi,2- At constant a and assum- 
ing the fraction of matter in the dense phase 
and the particle densities do not depend on 
temperature we then have 



,d£ 

dt 



Cv dT ^ j.dsd/3 
p At dp dt 



(S2 

p dt dp dt 



+(i-xp)r 



^1) 

dpdsi rr^pds2 
- + XpT 



dt dp 



+T 



dt 



{Ss + Sc) 



dt dp 
(11) 



volume average in the mixed phase 

Cv = (1 - Xv)cv,l + XvCV,2 • (12) 

The latent heat absorbed by a particle 
crossing the boundary between two phases 
in equilibrium is the temperature times the 
difference in entropy per p article between the 



two p hases, q = T[s2 — si] (iLandau &: Lifshitz 



3); 



19801 ): we recover this in the second term on 



the right hand side of Eq. (jlip . If any term 
in Eq. ([TT]) is negative, heat is evolved by this 
term which then heats the star and adds to 
the luminosity L. In particular the latent heat 
is a heating term for increasing density when 
the entropy per baryon of the dense phase 
is less than that of the low density phase - 
a situation which will arise when considering 
superfluid quark matter. 



For future reference we identify the terms 

(s2 - Si) is 



in Eq. ^ as follows. 



identified as the latent heat, (1 — Xp)T^^ 
is identified as the hadronic bulk contribu- 
tion, XpT^^ is identified as the quark bulk 
contribution, and the last term T(d/dt)[(S's + 
Sc)/p] is identified as the surface and Coulomb 
contribution. Further we often refer to these 
terms collectively as additional entropy pro- 
duction (or release) beyond what would be 
expected at constant density. 

We discuss in the following section how 
to calculate the bulk entropy density of the 
hadronic and quark matter phases in the 
relativistic mean field theory framework of 
the CfgQ equation of state. For now let us 
just remark that the entropy per particle for 
an ideal relativistic dege nerate Fermi gas is 
(jLandau Lifshitd[l98ol l 



(37r2 



She 



3 1 

-Tp'a 



0.02 



MeV Vfm 



P 



where the heat capacity is again a weighed 



(13) 

Taking the transition to quark matter as a 
transition between such gases - note that the 



7 



bag constant does not contribute to the en- 
tropy - and remembering that each hadron 
contains three quarks which further have a 
color degeneracy of three, the latent heat in 
Eq. (jlip is of the order of 



PQ 
fm~' 



fm~' 



daq/dt 
1057/107 yr 



(14) 



erg s 



where Tg = T/10^ K, and the rate of 
baryons making the transition to quark mat- 
ter, daq/dt, was (arbitrarily) scaled to an 
entire star being converted steadily over 10^ 
years. Unless the star is very hot or changing 
structure fast this is a very modest contribu- 
tion, and we further note that it is positive 
for a contracting star and so acts to cool the 
star down. However the sign is subject to the 
very rough assumption that both gases may 
be treated as ideal Fermi gases for the purpose 
of calculating their entropy, and it will change 
in a more detailed treatment. Specifically the 
quark phase may be color-superconducting 
with energy gaps as high as 100 MeV and 
corresponding critical temperatures of the or- 
der of 10^^ K. Below the critical temperature 
the quark specific heat and entropy density 
would be exponentially suppressed and could 
be ignored relative to the hadronic contribu- 
tion in Eq. (I14p which would then be a heating 
term. We shall return to both possibilities in 
later sections. 

We shall calculate the two bulk terms nu- 
merically in the following sections, so for now 
we just note that from Eq. ()13p the entropy 
per baryon decreases with increasing density. 
In a contracting star these terms will hence 
act as heating terms and locally be of the same 
order of magnitude as the release or absorp- 
tion of latent heat in Eq. (jl4p - but of course 
they contribute throughout the star and are 
therefore potentially far more important than 
the latent heat which is only significant in the 



mixed phase. 

Since the surface and Coulomb energies 
are related by £'s = 2Ec in equilibrium 
(iGlendennind boool ^. their contributions to 
the thermodynamic potential, and hence to 
the entropy, are similarly fixed in proportion, 
and we need only consider one of them here. 
Specifically the sur face part of the entrop' 



may be found from (jLandau fc Lifshitall98i 



Ss 



ons _ 

dT 



(15) 
(16) 



where r^s is the surface contribution to the 
thermodynamic potential 0,. The surface ten- 
sion for the quark-hadron interface, a, deter- 



phase ( 


Glendennins 


19921: 


Heiselbere et al. 


1993 


: Voskresenskv et al. 2003: 


Endo et al. 


2006 


). The surface tension of strangelets in 



vacuu m has been evaluated by lBerger &: Jaffe 



(|l987l ). 3ut the surface tension for the mixed 



phase remains essentially unknown. It is com- 
monly parameterized as being proportional to 
the difference in energy density between the 
two phases and the le ngth scale L ^ 1 fm of 
the strong interaction ( Glendenning 200d ) 



a 



K x[eQ- en] x L 



(17) 



Assuming for simplicity that K and L are con- 
stant we then get 



da 
df 



a 



cv,Q - cv,H 



eq - 
p-HSs + Sc) 



en 



3 es Cy,Q - Cy^H 
2 p eg - en 



(18) 
(19) 



The corresponding term in Eq. (jlip is then 
of the order of the surface energy per baryon 
times the ratio between thermal and total en- 
ergy density. It can therefore not be expected 
to contribute significantly, and this expecta- 
tion is confirmed by the numerical results. 
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3. Equation of State and Stellar Mod- 
els 

In our numerical work we have employed 
the rotating neutron star code developed by 
Weber and the G fsn equation of state used by 



Glendennin3 (fl992. 2000) and lGlendenning fc Weber] 
(|200ll ). Here we shall briefly describe each of 
these and the resulting stellar models. 

The G'^go equation of state treats the de- 
confined quark matter phase in a simple ver- 
sion of the bag model which ignores gluon in- 
teractions. The confined hadronic phase is de- 
scribed in terms of the mean field solution to a 
covariant Lagrangian that involves the baryon 
octet interacting through scalar, vector and 
vector-isovector mesons. The precise values of 
the coupling constants for the hadronic part 
of the mo del correspond to the first set of Ta- 
ble 5.5 of iGlendenninel (l200d 'l or Table 1 of 
Glendenning &: Webed (120011 ). For details we 
refer to these works where the resulting equa- 
tion of state as well as the underlying theory 
are carefully described at zero temperature. 
Finite temperature expressions for pressure, 
energy and particle densities can be found by 
reinserting the Fermi distribution in the phase 
space integrals of their zero temperature ex- 
pressions which we s hall not write explic itly 
here. We refer to e.g. Glendenning ( 199d ) for 
the full temperature dependent expressions 
(this reference also includes gluon interactions 
to first order which we ignore here). 
The entropy is calculated as 



1 

7t 



(20) 



contributions from thermally exited particle- 
antiparticle pairs. 

Figs. [T] and [5] show the resulting chemical 
composition and entropy. In the quark phase 
we note that u-quarks are suppressed initially 
giving the phase a negative net charge, and 
that as expected the entropy per baryon - 
though not per quark - is higher in the quark 
matter phase in the absence of any pairing 
phenomena. We have checked numerically 
that in the absence of pairing the entropy sim- 
ply scales linearly with temperature within 
the temperature range we shall need. 

The sum of the surface and Coulomb ener- 
gies has a maximum as a function of density 
around p = 0.5 — 0.6 fm~^ with a correspond- 
ing minimum in the related entropy in Fig. [3j 
The surface and Coulomb term in Eq. (jlip 
can therefore be either positive or negative 
and either heat or cool the star accordingly. 
The surface and Coulomb contribution to the 
entropy is negative with the total entropy re- 
maining positive, which confirms that a struc- 
tured phase has lower entropy than an un- 
structured one. 




0.001 



Total baryon density, p [fm 



using the finite temperature expressions out- 
lined above and using the highl y accurate 



Fig. 1. — Chemical composition of the GfgQ 
— -^qHation of state as a function of total baryon 
pubhcly available code described in| Miralles van RiQgji sity. Individual particle densities refer 
Um) to solve the Fermi integrals. In to the density in regions filled with the respec- 

tive phase. Xv is the (dimensionless) volume 
fraction of quark matter. 



keeping with our assumption that tempera- 
tures remain too low to significantly influ- 
ence the chemical composition we neglect 
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Hadronic phase combined 

Quark phase combined 

Protons 

Neutrons 

A particles 

u- quarks 

d- quarks 

s -quarks 




0.01 



Total baryon density, p [fm 



Fig. 2. — Entropy per baryon at T = 1 MeV 
in each phase, sh and sq, and per particle for 
each particle, Si, as functions of total baryon 
density for the GfgQ equation of state. Protons 
and A particles have high Si because they are 
so rare, but contribute little to sh for the same 
reason. 



-0.001 



+ -0.002 



-0.003 



-0.004 




Total baryon density, p [fm ] 

Fig. 3. — Sum of surface and Coulomb en- 
tropy as a function of total baryon density for 
the Gf^[] equation of state at T = 1 MeV 

We have calculated the structure of a se- 
quence of stars with rotation frequencies (as 
seen by an observer at infinity) between zero 
and the limiting mass-shedding Kepler fre- 
quency for the Gfgg equation of state. For 
th is purpose we use the perturbative method 
of lHartlel (|l967l l and lHartle fc Thornj (|l968l ) 
as implemented in the numerical code de- 



veloped by Weber which also solves self- 
consistently for the gener al relat i yistic Ke- 
pler frequency Qk - see IWebeii (119991 ') for 



the derivation of and (|\yeber et al.lll991 



Weber &: Glendennina 11992 ) for further de- 
tails. The sequence has constant total baryon 
number A = 1.87 x 10^^ and nonrotat- 
ing total gravitational mass M = 1.42 Mq. 
The Kepler frequency is then found to be 



= 6168 rad s 
of 1.02 ms. 



corresponding to a period 




7000 



ii [rad s 



Fig. 4. — Circumferential stellar radius in the 
equatorial and polar directions for rotating 
stars of total nonrotating gravitational mass 
M = 1.42 M0. 

Figs, m and [5] show a few properties of 
the models. The stars are significantly dis- 
torted by rotation and increase their equato- 
rial radius at the Kepler frequency by half the 
nonrotating radius (Fig. H]) while losing the 
pure quark matter core at O = 1400 rad s~^ 
(Fig. [5]). Since the stellar photon luminosity 
is proportional to the surface area and must 
match the energy flux emerging from the core, 
higher surface temperatures must also be ex- 
pected at low rotation frequencies for this rea- 
son alone. The distortion from spherical sym- 
metry depends on polar angle, however, and 
above the frequency at which the quark mat- 
ter core is lost the star actually contracts in 
the polar direction. 
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In Fig. Owe show the location of the phase 
boundaries between pure hadronic matter, the 
various geometries of the mixed phase and 
the pure quark matter phase. Here we use 
as the free variable the baryon number, a, 
contained within a surface on which the den- 
sity is uniform (i.e spatially but not tempo- 
rally constant). The Eulerian density change 
r can be positive or negative depend- 
ing on location while the Lagrangian deriva- 
tive (dp/di7)a is always negative, and so o 
is a more convenient variable for some pur- 
poses. In Fig. [5] borders are shown at the 
densities which correspond to each transition 
and a is scaled to the total baryon number, 
A. We note that a large fraction of the star 
is converted to quark matter as the star spins 
down. The pure quark matter core appears 
below Q = 1400 rad s~^ and grows to eventu- 
ally comprise 30 % of the stellar gravitational 
mass and 26 % of the stellar baryon number. 




Hadronic drops 

slabs 



0.4 0.6 
a/A 



Fig. 5. — Location of phase boundaries in 
stars of nonrotating mass M = 1.42 Mq. The 
locations are given as the fractions of the total 
stellar baryon number contained within a sur- 
face of spatially constant density correspond- 
ing to each transition. The inset shows the 
region where the phase boundaries almost - 
but not quite - cross. 



Within the framework of Hartk ( 19671 ) the 
changes in pressure and density with respect 
to a nonrotating star are second order effects 
in the rotation frequency , and as discussed 
by Weber &: Glendennine ( 19921 ) this remains 
true even at the limiting Kepler frequency, 
Qk, where the star would shed mass from the 
equator. The Lagrangian density change dur- 
ing spindown can then be reasonably approxi- 
mat ed by the order of magnitude estim ate (see 
also Fernandez &: Reisenegger (|200i)) 



fdp_ 



(21) 



Fig. [6] shows how the numerically calculated 
spin down compression rate compares to this 
estimate in the mixed phase - we have checked 
it at other positions as well. For stars with 
high spin frequencies and no pure quark mat- 
ter core Eq. (j2ip is a reasonable approxima- 
tion although an overestimate in some regions. 
It is off by approximately a factor of 0.1 below 
the frequency where the pure quark matter 
core forms. 




Fig. 6. — Frequency derivative of density 
scaled to the estimate in Eq. (|2ip at a/ A = 
0.33 as a function of spin frequency. The esti- 
mate is best at high spin frequencies but holds 
within a factor of 30 in all stars. 



Eq. ([2T]) can be used to give a rough es- 
timate of the effects of the additional en- 
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tropy production related to density change in 
Eq. (HH) 



(22) 



^ds _ Cydr ^ ^dsdp 
dt p dt dp dt 

p dt dp 



Taking the entropy per particle from Eq. ([13 
this implies that the bulk terms in Eq. ([TT 
are of the order of 

dt dp 3 Q| 



-10"^^r<? 



B 



1014 G 



(23) 



/O/6000 y 
V rad s~^ J 



fm 



erg s 



where we have taken $7k = 6000 rad s ^, as- 
sumed a standard dipole model for the spin- 
down so P = (5i9/3.2)2p-i with the spin 
period, P, measured in seconds and the mag- 
netic field, B in units of 10^^ G. As a rough 
estimate the total heating power from the 
bulk terms in Eq. (jlip in the absence of pair- 
ing and for a star of constant density and tem- 
perature with a baryon number of 10^'' is then 
of the order of 

W m^^T^( ^ AV ^/6000 y 



fm 



erg s 



(24) 



This includes only the bulk terms of 
Eq. ([TT]). but on general grounds one would 
expect this to be a reasonable approximation 
of the total additional entropy production. 
The latent heat has the same basic origin 
- the release of entropy with changing den- 
sity - and is therefore expected to be of the 
same order of magnitude as the bulk terms 
in Eq. (I24j) . In the following section we shall 
find that this is justified under most, but not 
all, circumstances. 



The estimated heating in Eq. (j24p is to 
be compared with the neutrino and photon 
luminosities which in the absence of pair- 
i ng phenornena c an roughly be estimated as 
(|Page et ahlbood l 



L., ~ lO^^To^ erg s"^ 



(25) 
(26) 
(27) 



where 'slow' refers to stars dominated by rel- 
atively inefficient neutrino emission processes 
such as the modified Urea cycle, bremsstrahlung 
or the pair breaking and formation process 
while 'fast' refers to stars dominated by the 
highly efficient direct Urea process. 

From Eq. (j24p we then see that while it may 
be possible to find combinations of tempera- 
ture, magnetic field and rotation frequency for 
which the additional entropy production dom- 
inates, such configurations may be difficult to 
realize in nature. For instance, in order for 
the heating term in Eq. (j24p to dominate the 
neutrino luminosities the star must be a hot 
millisecond magnetar. Such an object would 
be very short lived if it could be created in na- 
ture at all. For the heating term to dominate 
the photon luminosity - which is the more rel- 
evant term at late times and low temperatures 
- extremely high magnetic fields would again 
be required assuming a dipole model for the 
spin down. Such conditions - if they were to 
be realized - would be very short lived and 
therefore the more difficult to observe. 

From these simple estimates we then ex- 
pect that the release of entropy related to 
changing bulk density with spin down, the 
latent heat of phase transitions and changes 
in the surface and Coulomb energy of mixed 
phases (smaller still) will have little impact on 
the thermal evolution of compact stars and be 
difficult to observe. 
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4. Spherical Isothermal Cooling 

Estimates cannot replace detailed calcula- 
tions, and the ones discussed above further ig- 
nore all effects of superfluidity, the changing 
chemical composition and the rapid changes 
in structure at certain frequencies. We there- 
fore proceed in this section to explore these 
effects through numerical calculations taking 
the simplest possible approach to solve the 
thermal balance in Eq. ([3]). 

4.1. Numerical Setup 

Neutron stars are commonly taken to be 
isothermal after an initial thermal relaxation 
lasting 10-100 years in the sense that the red- 
shifted temperature T°° = Te* is constant 
below a thermally insulating layer in the outer 



with 



crust ( Gudmundsson et al. 19831 ) . Since we 



intend only to investigate the late thermal 
evolution of neutron stars, and as we do not 
expect the heating terms discussed here to 
change the isothermality by disturbing the 
thermal balance of any part of the star with its 
surroundings, we shall work within this same 
approximation. To show the full range of ef- 
fects, plots in the following also include young 
stars not expected to be isothermal and our 
results should only be taken as indicative at 
such early times. 

Assuming the star to be thermally relaxed. 



so the redshifted temperature T'^ 



IS 



constant below the outer crust, we then inte- 
grate Eq. (jll) to get a simplified equation for 
global thermal balance 



dt C ^ 



- L-) (28) 



C 



T OO 

^7 



daT— 
dT 

^ dsdp- 

da^e^* 
P 



(29) 
(30) 
(31) 



.(0) 



c 



47ri?V(r^"^)V(B)e2*^ (32) 



Tg - 0.001 5l4 



(34) 



where e^/p = Qu is the neutrino emissiv- 
ity, the surface temperature is given in units 
of 10^ K, Tsg = Ts/lO^ K and is 
surface temperature at zero magnetic field 
strength. We have used a fit appropriate 
for iron envelopes to the relation between 
the surface temperature, Ts, and the inner 
temperature below the heat-blanketing enve- 
lope, T, in which (714 i s the surface gravity 



in units of 10^'^ cm s ^ (Potekhin et al. 



1997 



Potekhin fc YakovlevI 120011 ). 

The heat conductivity in the crust, which 
derives mainly from the motion of electrons, 
becomes anisotropic in a stron g magnetic field 



(iPotekhin YakovlevI I2OOII : iGeppert et al 



20041 ) . It is slightly enhanced in the direction 
along the magnetic field lines and strongly 
suppressed in the direction orthogonal to the 
field lines. The heat blanketing relation and 
the surface temperature then becomes nonuni- 
form, TiocAB,e,g,T) = T^^\g,T)x{B,e,T). 
The photon luminosity must therefore be 
found by integrating the locally emerging fiux 
and it will depend on magnetic field strength, 
L^{B) = L^(0)F (B). For this purpose we use 
the fits given by lPotekhin Sz YakovlevI (I2OOII ) 
for xiB,9,T) and F{B) in a dipole magnetic 
field. These are good for magnetic fields below 
10^^ G, interior temperature between 10^ K 
and 10^ K and surface temperature above 
10^ K. F{B) reaches a minimum of ~ 0.7 
around B = 10^^ G and grows to about 2 at 



13 



B = 10^^ G for a star of internal tempera- 
ture T = 10^ K. The location of the minimum 
moves to lower magnetic field strength for 
lower internal temperature and the change in 
photon luminosity can delay or accelerate the 
cooling at late times accordingly. The figures 
in Sect. 14.21 show the effective surface tem- 
perature at infinity = F{B)'^/^e'^ . A 
strong magnetic field in the inner crust can 
also lead to anisot r opic heat flow there (see 
e.g. lOeppert et al" ( 20041 ) and Aguilera et al 



(|2008l X^. l)ut these effects cannot be included 
here. 

For given assumptions about the neutrino 
emissivity and superfluid gaps Eq. is 
solved along with the spin down model to give 
the surface temperature at infinity func- 
tion of time, rotation frequency and magnetic 
field strength. At each time step Eqs. ([29]) 
to are solved to update the stellar struc- 
ture at the appropriate rotation frequency. 

Our equations of thermal balance assume 
spherical symmetry. This effectively treats 
rotation as a perturbation whose only effects 
are to change the size and chemical compo- 
sition of the stars and to provide additional 
terms in the thermal balance. This approach 
ignores the effects of nonradial heat fiows 
and nonspherical perturbations of the metric 
which can be significant when the star is not 
spherically symmetric at very high spin fre- 
quencies (see Fig. H|). Two-dimensional cool- 
ing calculations are beyond our scope, but 
they have been pe rformed at constant ro - 
tati on frequency b y Schaab Sz Weigel ( 19981 ) 
and iMiralles et al.l (|l993l ) . For stars rotating 
at a large fraction of their Kepler frequency 
these authors found significant effects on the 
thermal evolution during the nonisothermal 
epoch, and polar temperatures up to 31% 
higher than the equatorial temperatures even 
for internally isothermal stars. The recent 



work of lCeppert et al. ( 2004 ): Aguilera et al 
;'urther includes the influence of large 
magnetic fields in the crust and Joule heat- 
ing by magnetic field decay in two dimensions. 
This results in nonradial heat flow and signif- 
icant heating from magnetic fleld decay. We 
focus on the effects of deconfinement during 
spin down however and shall not pursue these 
aspects here. 

Unless otherwise stated we shall use the 
equatorial radius in Eqs. ([32]) to ([M|) . We 
have performed calculations with the polar 
radius too and will show one such example. 
At high rotation frequency the surface tem- 
perature is then higher and responds less to 
changes in rotation frequency. We expect 
results in a two-dimensional code would be 
intermediate between these and such a calcu- 
lation would be a significant improvement on 
what is presented here. With this in mind we 
use the equatorial radius to show the most 
pronounced effects of spin down. However, 
we also note that the strongest effect we have 
found from spin down stems from the for- 
mation of a pure quark matter core. At the 
rotation frequency where this happens the po- 
lar and equatorial radii are very similar (see 
Fig. HI) and a spherical approximation is rea- 
sonable. 

Since the stellar moment of inertia / and 
radius also change with rotation frequency 
we modify the spin down model described 
in the previous section to include such ef- 
fects. The spin do wn is then determined by 
(iGlendennindbood ^ 



K 
T 



1 + 



I' 



dl 



K = {2/3)R^B'^ sin^a 



(35) 
(36) 



where a is the inclination angle between the 
magnetic axis and the rotation axis and m — 1 
is the multipolarity - usually m = 3 for mag- 
netic dipole braking or m = 5 for gravitational 
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quadrupole radiation. We take sin^ a = 1 and 



use the canonical value m 
spin down. 



3 for a dipole 



As for the neutrino emissivities and super- 
fluid properties of the neutron star we are in- 
terested in the most important effects only 
and aim for transparency in the model. In the 
quark phase we include the direct and modi- 
fied Urea cycles as wel l as bremsstrahlung us - 



ine the emissivities in Blaschke et al. I Sw 



The US-branch of the direct Urea cycle fulfills 
the triangle inequality, 



(37) 



where kp- is the Fermi momentum of the 
strange-quark, up-quark and electron respec- 
tively. It is thus allowed throughout the quark 
phase, but its ud-counterpart is forbidden as 
it cannot conserve both energy and momen- 
tum. This difference stems from the assumed 
strange quark mass of 150 MeV. 

In the hadronic phase we again include the 
direct and modified Urea cycles as well as 
br emsstrahlung u s ing th e emissivities given 



m 



Yakovlev et al 



((20011). We ignore the ef- 
fects of hyperons here because of their low 
abundance. In neutron stars the direct Urea 
reaction is often allowed only at very high 
densities because it cannot fulfill both energy 
and momentum conservation unless the pro- 
ton fraction is above the value where both 
charge neutrality and the triangle ine quality 
can be observed ( Lattimer et al. 199ll ). It is 
critical to note that this is not so in hybrid 
stars with a mixed phase. The mixed phase 
is possible precisely because charge neutrality 
does not have to be observed locally, and it is 
favorable because the hadronic phase lowers 
its nuclear symmetry energy by increasing the 
proton fraction (Fig. [1]). Hence the hadronic 
direct Urea cycle is active in the mixed phase, 
and as it is about three orders of magnitude 
faster than the us-branch of the quark direct 



Urea cycle it controls the cooling of hybrid 
stars with a mixed phase unless reduced by 
pairing effects. 

Additional neutrino processes should be 
included in a more sophisticated treatment - 
particularly in the crust - but as stated above 
we are mainly interested in the general prop- 
erties of the additional entropy production 
and the changing chemical composition and 
we shall here leave out such terms. 

If there is any attractive interaction among 
particles in a degenerate Fermi system they 
will pair, and the resulting superfluidity has 
important consequences for the thermal prop- 
erties of neutron star matter and the thermal 
evolution of neutron stars. Pairing in general 
delays the cooling because it suppresses most 
of the neutrino emissivities and enhances the 
heat capacity at temperatures just below the 
critical. However it also opens additional neu- 
trino emission processes, suppresses the heat 
capacity far below the critical temperature 
and enhances the thermal conductivity, and 
it may therefore also ac celerate the coo l ing a t 
certain epochs (see e.g. Yakovlev et al. ( 19991 ) 
for a detailed account). 

There is considerable uncertainty concern- 
ing the rel evant superfluid ph ase in quark 
matter - see Alford et al. ( 20081 ) for a compre- 



hensive review of the effects of pairing among 
quarks and the possible range of phases. For 
this first investigation of consequences for the 
thermal evolution of a quark-hadron phase 
transition during spin down we assume a sim- 
plified but physically transparent model for 
quark matter pairing which essentially corre- 
sponds to_2ainn^in__the color- fiavor l ocked 
phase dAlford et al.l [l999l : lAlfordI l200lh . In 
this model each fiavor participates equally, 
the gap Aq_o is independent of density and 
the critical temperature is Tc = 0.72Aq^o (see 



Schmitt et aD (|2002l )). We shall consider 



wide range in Aq^q partly because it is in- 



15 



teresting for the present calculation, partly 
because very low gaps may be realized for 
high strange quark mass as this implies un- 
equal quark Fermi momenta in the unpaired 
phase. The color-flavor locked phase is gener- 
ally understood to be so strongly suppressed 
in all its thermal properties as to be virtually 
inert with respect to cooling. Since we shall 
also consider small gaps and since the latent 
heat derived from the phase transition de- 
pends strongly on pairing in the quark phase 
this will not always apply here, but for strong 
quark pairing the thermal evolution is gen- 
erally controlled by the hadr onic phase. We 
furthe r note that as sho wn by Jaikum ar et al.l 
(|2002l ) (and stressed by lAlford et al . (2008)) 
Goldstone modes and their related neutrino 
emissivities and heat capacities are not expo- 
nentially suppressed, scaling instead as T^^ 
and respectively. They will therefore dom- 
inate in color-flavor locked quark matter at 
low temperatures, but they are exceedingly 
small and do not influence our results. 

The effects of superfluidity relevant for our 
purpose are that below the critical tempera- 
tures it suppresses the neutrino emissivities 
and entropies of the participating particles 
while also allowing additional neutrino emis- 
sion through Cooper pair breaking and forma- 
tion. The speciflc heat is enhanced just below 
the critical temperature and suppressed at 
lower temperatures. These effects are incor- 
porated through control functions, i?(T, Tc), 
which multiply the relevant quantities. They 
depend only on the pairing channel, temper- 
ature and gap size, and may therefore be cal- 
culated given expressions for the gap size mo- 
mentum dependence. 

In the quark matter phase we use the 
simplest possible control functions and sup- 
press the direct Urea emissivity by a fac- 
tor of e"^'^/'^ and the modifled Urea and 



bremsstrahlung emissivit y by a fact or of 

g-2AQ/T^ whe re (followingl &teinereLalJ (jiool)) 
Aq = Aq,o\/1 - {T/TcY is the pairing gap 
at the local temperature, T = T°°e~^ . 

The specific heat and entropy of super fluid 
particles are also modified at temperatures be- 
low the critical temperature. In the quark 
matter phase we use the fit to the ^Sq con- 
trol function f o r the specific heat given by 
Yakovley et al.l (119991) . We further fit the re- 
sults of iMiihlschlegell (|l959l ^ for the entropy 
suppression to obtain 

Ssf = 0.95so X (T/Tje-^^/^ 



0.43 + 3.82(r/rc) - 1.41 (T/Tc 



(38) 



Where sq is the nonsuperfluid entropy. This 
expression is applied to reduce the entropy 
of both quarks and hadrons (thus here ne- 
glecting the difference between different pair- 
ing states of neutrons). The entropy dif- 
ference between the superfluid and the nor- 
mal phase depends on temperature as As oc 
— (1 — T/Tn) close to the critical temperature 
(jLandau &: Lifshitd [l98ol l with So becoming 
negligible below 0.2Tc. This has the effect of 
gradually turning the latent heat in Eq. (jlip 
from a cooling term into a heating term, but 
it also suppresses the bulk terms in Eq. ([TT]) . 

In the absence of detailed calculations of 
control functions for the entropy of superfluid 
quark matter phases in the literature we have 
employed Eq. (j38p in the quark matter phase 
though it is strictly relevant only for BCS su- 
perfluidity. Eq. ([38]) is, however, consistent 
with the exponential suppression of entropy 
below the critical temperature which would 
be expected in any such calculation, and we 
do not expect that our results are sensitive to 
this choice. 

Pairing further opens the important pos- 
sibility of neutrino emission by Cooper pair 
breaking and formation in the quark phase 
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and we incl ude the associated emissivit y as 
described in Jaikumar &: Prakash ( 200ll ) us- 



ing the fit to the ^ Sn control function in 



Yakovlev et "aD (|l999l l. 



Among the nucleons pairing is predicted for 
neutrons in the ^Sq and channels and for 
protons in the ^Sq channel. The gaps can be 
calculated self-consistently but results are un- 
certain and vary greatly both in terms of the 
maximum gap size and in terms of density 



depen d ence - see e .g. Lombardo Sz Schulzj 



(|200lh : IPage et al.l (|2004l . bood l for collec- 



tions of these results. We shall here use 
gaps with a phenomenolog ical momentum 



depen d ence as suggested by iKaminker et al 



(I2OOII'). lAndersson et al.l (|2005l ) and recently 
Aguilera et al.l toO^ ) 



Aik 



An 



oj 



{kF,N - ko, 
{kF,N 



k2'' 



{kF,N - k2f + kl 



(39) 
(40) 



where kp^N is the Fermi momentum of the rel- 
evant nucleon, N = n,p. Our choices for the 
paramet ers An and kj corresp ond to sets a, e 
and h of Aguilera et al. ( 20081 ) (see that work 
for references to the model calculation under- 
lying these fits). Eq. ()40p is valid only in the 
range ko < k < k2 with A{kF,N) = outside 
this range, and where the gaps for ^Sq and 
^P2 pairing of neutrons are both nonzero we 
use the largest of the two gaps. We further 
use the relations between cr itical tempera- 



ture a nd pairing gap listed by I Yakovlev et al 
(|l999l ). 



In the hadronic phase neutrons and protons 
may pair simultaneously in different chan- 
nels and so the resulting calculations for the 
control functions can be quite involved. We 
use the fits to numerically calculated control 
functions for each process and ea c h pai ring 
channel compiled in jyakovlev etHI JigS) (or 



Yakovlev et al.l ( 200ll ) - see these works for 
detailed references). Where protons and neu- 
trons pair simultaneously we use the com- 
bined factors listed in these works if available 
and if not then the smallest of the two inde- 
pendent control functions. 

The hadronic pair breaking and forma- 
tion process can be very powerful and may 
dominate the thermal evolution when pair- 
ing first sets in, bu t it wa s recently shown 
by iLeinson &: Perez (|200fil ^ to be strongly 



suppressed in the singlet state channel by 
approximately a factor of 10~^ relative to 
the ^P? channel ( see al so the recent work of 
Steiner Sz Reddv ( 20091 )). In the hadronic 
phase we therefore include pair breaking and 
formation only for neutrons at high densities 
where they pair in the '^P2 channel and for 
this we again us e the emissivity and contro l 
function given in lYakovlev et all (|l999l . [200ll ). 



4.2. Numerical Results 

Figures [7] to [l2] explore consequences for 
the thermal history of hybrid stars of includ- 
ing time dependent structure and entropy pro- 
duction by spin down compression in the cool- 
ing calculations with the approximations dis- 
cussed above. They show the surface temper- 
ature and additional entropy production at in- 
finity as functions of time, magnetic field and 
quark pairing gap energy. 

In Fig. [7] we show the thermal evolution of 
stars with quark pairing gap Aq^q = 10 MeV, 
a range of magnetic fields and initial spin pe- 
riod around 1 ms. The extremely rapid initial 
rotation is interesting because it means the 
stars start out with no pure quark matter 
core and undergo drastic changes in structure 
as a pure quark matter core develops around 
Q = 1400 rad s"^ We see this in Fig. [7| as 
a short period of increasing temperature at 
the time corresponding to this angular veloc- 
ity for a given magnetic field. During this 
short period the heating term from the addi- 
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Fig. 7. — The effect of spin down compression 
on cooling curves for stars with initial spin 
frequency restart = 6000 rad s~^ and constant 
magnetic field as labeled for each curve. The 
quark core is assumed superconducting with 
^(9,0 = 10 MeV. The increasing quark matter 
fraction produces slower cooling and a jump in 
temperature with the appearance of the pure 
quark matter phase. 

tional entropy production actually dominates 
which is shown in more detail in Fig. [10] dis- 
cussed below. Further, the changing neutrino 
emissivity and heat capacity with the intro- 
duction of more quarks in the star and the 
reduction in stellar surface the pure 

quark matter core develops all affect the ther- 
mal evolution. This implies somewhat slower 
cooling - with our specific choice of param- 
eters - and higher surface temperature. For 
the internal temperature, T, we have found 
essentially a transition between two otherwise 
similar cooling tracks; one for stars with no 
pure quark matter core and one for stars with 
a fully developed quark matter core. The ef- 
fective surface temperature also depends on 
magnetic field strength through F(B) how- 
ever. For this reason such an easy interpreta- 
tion is difficult from Fig. [7] alone. 

The magnetic field strength determines the 
spin down and affects the heat blanketing re- 



Fig. 8. — Variation of surface temperature at 
infinity with magnetic field for stars at ages 
10^, 10^ and 10^ years from above. Thick 
lines have initial spin frequency Ostart = 
6000 rad s~^ while thin dashed lines start 
with restart = 600 rad s~^. Thin continuous 
lines use the polar rather than the equato- 
rial radius. The quark core is assumed su- 
perconducting with Aqa = 10 MeV. We note 
a slow increase in temperature for old stars 
and a sudden jump at the magnetic field cor- 
responding to significant spin down at specific 
ages. The suppression in temperature at low 
magnetic field and increase at high field is due 
to the effects of the magnetic field of the heat 
blanketing relation. 

lation. Fig. [8] shows how the temperature 
varies at specific ages with the magnetic field 
strength (which is kept constant in time it- 
self). This is shown for initial spin frequen- 
cies of 6000 rad s~^ and 600 rad s~^ and using 
the polar radius in addition to the equatorial 
when calculating the surface temperature. 

For low initial spin frequency we find no 
significant effects of spin down in Fig. [8l The 
change in temperature at a specific age with 
increasing magnetic field is due to the effects 
of the magnetic field on the heat blanketing 
relation and photon luminosity. If these ef- 
fects were left out the curves for low initial 
spin frequency would be almost constant. 
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If the initial rotation frequency is suffi- 
ciently high to exclude the pure quark mat- 
ter phase it is a different matter. A strong 
magnetic field initially suppresses the surface 
temperature of young stars with high rotation 
frequencies, but it then gives a sharp jump 
at the magnetic field strength which ensures 
that a pure quark matter core is formed at a 
particular age. For older stars there is also 
a slight increase in surface temperature be- 
fore the formation of a pure quark matter core 
when the equatorial radius is used. If the ef- 
fects of the magnetic field on the heat blan- 
keting relation were left out the curves at high 
initial spin frequency become nearly constant 
for B > 10^^ G and for fields too weak for 
spin down to set in at the particular age. In 
this case there would still be a jump in tem- 
perature at the magnetic field strength corre- 
sponding to formation of a pure quark matter 
core for a specific age however. 

Replacing the equatorial radius with the 
polar in Eqs. ([52]) and (thin continuous 
lines in Fig. [8]) we find higher surface temper- 
atures before the introduction of pure quark 
matter in the core and still a clear jump in 
temperature after. The gradual increase in 
temperature before the jump found in the 
other curves is absent except for old stars. 
This can be understood if we remember that 
the polar radius is smaller than the equato- 
rial and less sensitive to spin down at high 
rotation frequencies. We have checked that a 
similar pattern may be seen in plots of tem- 
perature versus time using the polar radius. 
We expect that full two-dimensional calcula- 
tions would give results which are at high ro- 
tation frequencies intermediate between what 
we have found using a spherical approxima- 
tion. The similarity between curves with 
polar and equatorial radius at the jump in 
temperature shows that our calculations are 
insensitive to the approximation at the spin 
frequency where the pure quark matter core 



forms and where we find the strongest effect 
of spin down. 
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Fig. 9. — Numerical value of the entropy pro- 
duction as defined in Eq. (|30p as a function of 
time for a star with magnetic field B = 10^^ G 
and quark pairing gap Aq^q = 10 MeV. Wrot 
is split into contributions with different phys- 
ical origins as explained in the text. The sur- 
face and Coulomb term changes sign and the 
negative of this term is shown as well. The 
latent heat and bulk terms are dominant, and 
the spike is caused by the appearance of a pure 
quark matter phase. 

In Figs. [9] and [10] we explore the importance 
of the heating term in a little more detail. 
Fig. [9] shows the numerical value of W°° and 
its various contributions as a function of time 
for a star with pairing gap Ag^o = 10 MeV 
and magnetic field B = 10^^ G. The total en- 
tropy production, 1^°°, is split into contribu- 
tions corresponding to (integrals of) the terms 
in Eq. pT|) . VF^^g^^ corresponds to the latent 
heat, W^H^dronic bulk ^° hadronic bulk term, 
^Qxiark bulk *° quark bulk term and W^q 
to the surface and Coulomb term of Eq. (jlip . 
Leptons do not participate in the deconfine- 
ment transition, so their entropy is assumed 
continuous across the phase boundary. They 
are separated from the other particles in the 
form of T^L^ton m Fig. B 
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Fig. 10. — The total entropy production 
scaled to the combined neutrino and photon 
luminosity for stars with quark paring gap 
^Q,o = 10 MeV and magnetic fields as in- 
dicated just above each curve. The different 
line styles are meant only to help guide the 
eye. Note that at the peaks reaching above 
1 the entropy production actually dominates 
the cooling terms and the temperature in- 
creases. Also note that this is only possible 
for stars changing structure rapidly with the 
introduction of a pure quark matter core. 

The bulk and lepton contributions always 
act as heating terms while the surface and 
Coulomb term changes sign and turns into a 
heating term with the onset of hadronic su- 
perfluidity around age 4000 years with the 
hadronic bulk term simultaneously falling 
away. The quark bulk term is very small 
because for a pairing gap of 10 MeV the 
quark matter is strongly superconducting at 
all times except the very early and the quark 
entropy is therefore suppressed. This also 
implies, however, that hadrons making the 
transition across the phase boundary essen- 
tially release all their entropy, and so the la- 
tent heat, which changes sign and becomes a 
heating term with the onset of quark superflu- 
idity, is therefore quite significant and actu- 
ally comes to dominate around the time when 
the hadronic bulk term falls away. The lep- 



ton term and the surface and Coulomb term 
are relatively small but significant when the 
bulk terms are eliminated through the onset 
of pairing. The peak in the total around age 
20000 years corresponds to the appearance of 
a pure quark matter phase in the core and 
the ensuing rapid changes in structure, and 
it contributes to produce the short interval of 
increasing temperature seen in Fig. [71 

Fig. do] shows the total entropy production 
scaled to the total neutrino and photon lumi- 
nosity, W°°/{L^ + -L^), for the same set of 
parameters which gave the cooling curves in 
Fig. [71 The strong temperature dependence of 
the neutrino and photon luminosities ensures 
that the scaled heating term does not increase 
above unity where heating and cooling exactly 
balance - except at the peaks caused by the 
appearance of a pure quark matter core. We 
also note that as in Fig. [71 the magnetic field 
must be above ~ 10^ G for to have any 
discernible effects. 
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Fig. 11. — The total entropy production 
scaled to the combined neutrino and pho- 
ton luminosity for stars with quark paring 
gap as indicated in the legend and restart = 
6000 rad s~^. The magnetic field is constant 
at 10^^ G. The entropy production changes 
sign and becomes negative for low pairing 
gaps and thus acts powerful cooling term 
during some epochs. 
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While the latent heat for a transition to 
a strongly superfluid core is a heating term, 
the transition to a nonsuperfluid core with 
higher entropy per baryon cools the star. This 
is ihustrated in Fig. [11] where we plot the 
total heating term scaled to the combined 
neutrino and photon luminosities at constant 
magnetic field for a range of quark pairing 
gaps Aq^o- Here we see how the entropy pro- 
duction changes from a heating term to a cool- 
ing term and may be quite dominant at cer- 
tain times depending on the choice of stel- 
lar parameters. This can also be observed 
in Fig. [12] where we plot surface tempera- 
ture versus age with constant magnetic field 
for a range of quark pairing gaps. The on- 
set of superfluidity initially delays the cooling 
and then accelerates it. When the pure quark 
matter core forms there is a drop in tempera- 
ture for very low pairing gaps and a jump for 
high pairing gaps. In fact, however, the ma- 
jor differences between the curves in Fig. [12] 
are due to the differences in the thermal evo- 
lution of superfluid and nonsuperfluid quark 
matter rather than spin down related effects. 
Old stars with a nonsuperfluid quark matter 
core cool more slowly than those without a 
quark matter core. 

A rich picture has now emerged from the 
coupling between the thermal evolution and 
spin down compression with deconflnement. 
The spin down may change the thermal evo- 
lution of hybrid stars by heating and cooling 
them at various ages and by changing the 
structure and chemical composition - or it 
may be entirely inconsequential depending on 
quark superfluidity and initial spin frequency. 
The effects we found turned out to depend 
strongly on the inclusion of quark pairing and 
the timing of the appearance of the pure quark 
matter phase in the stellar core. These are 
subject to the specific assumptions made con- 
cerning the equation of state, stellar baryon 
number and pairing regime, so it is essential to 
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Fig. 12. — The effect of spin down com- 
pression on cooling curves for stars with ini- 
tial spin frequency Ogtart = 6000 rad s~^ and 
quark pairing gap Aq^q as labelled for each 
curve. The magnetic field is constant at 
10^1 G 

stress that the results discussed in the present 
section illustrate the general point that spin 
down and thermal cooling may be interde- 
pendent, rather than providing quantitatively 
reliable predictions. 

Bearing this in mind we show in Fig. [13] 
how calculations for selected pairing gaps and 
magnetic field strengths measure up to obser- 
vational data on the thermal state of isolated 
neutron stars (kindly supplied by Dany Page) . 
This figure is shown partly as a consistency 
check for our calculations and partly to com- 
pare the size of the effects we have found with 
the accuracy of actual measurements. We 
stress that the observational data are consis- 



tent with static cooling models (|Page et al 



20041 ) and that we do not suggest that the 
effects of spin down are necessary to explain 
them. Further the effects we have found oc- 
cur after thermal relaxation only for mag- 
netic field strength below ^ 10^^ G. For those 
sources in Fig. [13] for which the field is k nown 
it is above that value ( Pons et al. 200?! ). 
Fig. [13] explores a range of quark pairing 
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gaps which in our pairing scheme are den- 
sity independent. In nature the pairing gap 
would not span a range this broad, but it 
could be density dependent and variations be- 
tween stars in mass and density would then in- 
troduce variations in their thermal evolution. 
Similarly stars with different masses would 
have a different phase structure and acquire 
pure quark matter cores at different rotational 
frequencies thus further complicating the pic- 
ture. 

Because the hadronic direct Urea mecha- 
nism and the neutron pair breaking and for- 
mation mechanism are active in the mixed 
phase our models are consistent only with 
the cold sources and the effects of spin down 
do not change this conclusion. That these 
mechanisms are the cause of the generally low 
temperatures in our models can be seen if 
they are artificially left out. The two dot- 
ted lines in Fig. [13] show - for illustration 
only ~ that temperatures are much higher if 
these mechanisms are suppressed. Such a sup- 
pression could occur for certain quark pair- 
ing schemes where pairing forces the quark 
phase to become electrically neutral or pos- 
itive rather than negative, thereby reducing 
the proton content in the hadronic phase be- 
low the threshold for direct Urea. A self- 
consistent inclusion of such effects are beyond 
the scope of the present investigation. We 
stress again that the inconsistency with hot 
sources seen in this figure is subject to the 
choice of stellar and thermal models made for 
the specific purpose of studying the effects 
of deconfinement during spin down. Thus 
it should not be seen as indicating a general 
breakdown in cooling theory. 

The true signal of deconfinement we have 
found is in the transition between cooling 
curves with and without pure quark matter 
cores, and in the brief period of increasing or 
rapidly decreasing temperatures thus caused 
for certain stellar parameters. As can be seen 



in Fig. [13] the change in surface temperature 
caused by deconfinement is smaller than the 
error bars on the observational data and it will 
be difficult to test observationally on the basis 
of data relating only temperature and time. 
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(2) IE 1207-52 

(3) RX ja002-»6246 

(4) PSR 0538+281 7 

(5) PSR 1055-52 

(6) RXJ 0728.4-3125 

(7) PSR 0633-45 

(8) PSR 1706-44 

(9) PSR 0656+14 
CIO) PSR 0633+1748 
(11) RXJ 1856-3754 




Fig. 13. — Comparison between observa- 
tional data and cooling calculations with spin 
down. The stars have initial spin frequency 
^^start = 6000 rad s~^, magnetic field strength 
B = 10^^ G and quark pairing gaps as indi- 
ca ted. The observa tional data can be found 
m IPage et al.l (j2004l ) and sources are identified 
by numbers right above or below their temper- 
ature error bars. The two dotted lines illus- 
trate slow cooling with artificial suppression 
of neutrino emission by direct Urea and neu- 
tron pair breaking and formation. They have 
Aq,o = 10 MeV (thick hne) and 0.01 MeV 
(thin line). Including these mechanisms only 
cold sources can be explained. The jump in 
temperature at the formation of a pure quark 
matter core is smaller than the error bars on 
the observational data. 

5. Discussion 

Our intention with the present work was 
to explore a possible connection between the 
spin down and thermal cooling of hybrid stars 
with a deconfined quark matter core - the 
appearance of which might be expected to in- 
fluence theoretical cooling tracks. Specifically 
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we were interested in the latent heat of the 
phase transition and the drastic changes in 
structure and chemical composition resulting 
from the increasing density with spin down. 
The general formalism worked out in Sect. [2] 
would also be relevant to models with no 
phase transition, however, as the bulk terms 
in the spin down derived entropy production 
might also be important in such cases. While 
estimates of the importance of the spin down 
derived entropy production did not give cause 
to expect strong signals of deconfinement in 
the thermal evolution, numerical calculations 
revealed clear effects of the changes caused by 
the appearance of a pure quark matter core. 
This signature is related to changes in radius, 
chemical composition, structure and under 
certain circumstances the additional entropy 
production both in bulk and in the form of 
latent heat. 

Our numerical work was carried out using 
a sophisticated and reliable code with respect 
to the stellar structure solving Hartles pertur- 
bative equations for the structure of rotating 
compact stars self-consistently. The thermal 
evolution, however, was treated in a spherical 
isothermal approximation with a somewhat 
ad hoc approach to the treatment of super- 
fluid pairing in the quark matter phase. Still 
we believe our work is sufficiently detailed 
to demonstrate the general point that impor- 
tant signals may be derived from the inter- 
play between spin down and cooling of com- 
pact stars, which could in the long run help 
answer pressing questions about the state of 
matter at high densities. The signals we have 
found do not dominate the general thermal 
evolution of hybrid stars but they do com- 
plement the standard picture. By correlating 
temperature with magnetic field strength and 
spin frequency they may also help break the 
degeneracy between models relating only tem- 
perature and time. 



The signature of deconfinement found here 
is below the present observational sensitiv- 
ity and not of sufficient strength to set apart 
the cooling curves with temperature versus 
time for hybrid stars. The correlation between 
temperature and magnetic field strength pro- 
vides a possible alternative. To test such a 
correlation it would be necessary to obtain ac- 
curate temperatures for a number of stars of 
approximately the same age for which the spin 
down could also be detected by the emission of 
pulses in either radio or X-ray. It might then 
be possible to test for features in the relation 
between temperature and magnetic field - or 
equally interesting, the rotation frequency it- 
self - which as demonstrated here could result 
from a strong phase transition if the initial ro- 
tation frequency was sufficiently high to have 
spun out the high density phase. 

Contemplating these prospects it is im- 
portant to consider that alternative effects 
not treated here might determine the rela- 
tion between spin down and thermal evolu- 
tion. Most pressingly the h eat from magnetic 
field decay was employed bv lPons et al.l (j2007l ) 
to explain the observed correlation between 
temperature and magnetic field and it may 
well drown out any other effects ~ although 
as noted by the authors the magnetic field 
decay itself has not yet been independently 



demon stra ted. Recent work bvlAguilera et al 



(|2008l l and lAguilera et al.l ([200^) on the cool- 
ing of magnetized neutron stars in two dimen- 
sions also highlighted the importance of the 
magnetic field strength, geometry and possi- 
ble decay for the thermal evolution of neutron 
stars. These authors find strongly anisotropic 
surface temperature distributions and possi- 
bly an inverted temperature distribution with 
hot equatorial regions for middle aged stars. 
They further show that the effects of magnetic 
field decay and Joule heating can dominate 
the thermal evolution at strong and interme- 



23 



diate field strengths - even to the point that 
the effects of the direct Urea process may be 
hidden by this heating term in magnetars. 



and it may well be lower than discussed here 
thus changing the timing of spin down related 
effects. The strongest signature of deconfine- 



Th e work of lReiseneggeij (|l995l ) and [Fernandez &: 
is also most important. These authors 
consider the second term on the right hand 
side of Eq. ^ which is shown to give rise to 
so-called roto-chemical heating as the weak 
interactions required to maintain chemical 
equilibrium are un able to keep pace with the 
increasing density. [Fernandez &: Reisenegger 
Found that old millisecond pulsars 
reach a quasi-equilibrium in which the photon 
luminosity is determined entirely by the spin 
down power and remains much higher than 
otherwise attainable. Their work considered 
only nucleonic equations of state but similar 
results should clearly be expected for hybrid 
stars just as part of the effects demonstrated 
here should be expected to show up in nucle- 
onic stars. 

Given the possible importance of such al- 
ternatives and the shortcomings of the present 
study discussed above, a more complete treat- 
ment of the interplay between the spin down 
and thermal evolution of neutron stars seems 
desirable before strong assertions can be 
made concerning the specific shape of cool- 
ing tracks. As well as including all possible 
energy sources and sinks this should be ex- 
tended to consider a wider range of stellar 
masses, equations of state, phase transitions 
and pairing regimes. A two dimensional code 
could further treat the influence of the mag- 
netic field on heat flows in the inner crust 
and the effects of nonspherical heat flows for 
rotationally perturbed stars - aspects which 
are all independently well described in the 
literature. 

It should further be noted that the rota- 
tion frequency at which the pure quark mat- 
ter core appears depends strongly on the stel- 
lar baryon number and the equation of state. 



gjtgjggh^gj^ i here depends entirely on the for- 



mation of a pure quark matter core where 
none existed previously. If, on the other hand, 
a pure quark matter phase is present in the 
core of stars at even the highest rotation fre- 
quencies - as is the case for other equations of 
state - we would not expect equally strong sig- 
nals from deconfinement with spin down com- 
pression. 

It is possible that circumstances not con- 
sidered here could provide a more pronounced 
link between rotation and cooling. If for in- 
stance the magneti c field va r ies oy er time - 
and as discussed by Geppert ( 20061 ) this may 
well be the case - the timing of any strength- 
ening or decay of the magnetic field may 
provide entirely different spin down histories. 
Alternatively the temperature of stars being 
spun up in accreting binaries should to some 
extent be determined by the changing chemi- 
cal composition and this could provide an al- 
ternative approach. 

Additional energy release would also be 
expected if the phase transition cannot main- 
tain thermodynamic or hydrostatic equilib- 
rium as was assumed here. As the most ex- 
treme example of this, we have entirely ig- 
nored the instabilities and corequakes which 
might accompany the appearance of a pure 
quark matter phase in a s tar far from equilib- 
rium (|Zdunik et al.l I2OO6I I . Such events could 
release large amounts of energy and signif- 
icantly change the thermal evolution by re- 
heating the star. 

While many essential issues thus remain 
poorly explored we hope to have demon- 
strated the possible usefulness of considering 
the connection between the spin down and 
thermal evolution of neutron stars. Although 
the results of such work may not be imme- 
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diately applicable to observational tests, we 
believe they could prove most valuable in the 
long run. 
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A. Appendix A, Derivation of the Energy Balance in Mixed Phases 

To derive Eq. (|3D from Eq. we proceed analogously to Thorne ( 19661 ) and Weber ( 19991 ) but 
work in terms of baryon number a instead of radial coordinate. We first define the nuclear energy 
generation rate per baryon as the rate, q, at which rest mass is converted to internal energy as 
measured by an observer in the shell and at rest with respect to the shells reference frame 



d?Tif 



d(proper time) 



dniB 
'e*dt 



(Al) 



where ttib is the average rest mass per baryon and we work in units with. The first term on the 
right hand side of Eq. ^ is the amount of rest mass converted to internal energy in the entire shell 
which is then q6ae^dt. 

Next we note that by expressing the shells volume in terms of its particle density, V = {V/ 5a)5a = 
p/5a, the second term in Eq. ([1]) can be written 



dW = -Pd {^6a] + adS + d^c 



(A2) 



The third term represents the difference between the rates at which energy enters and leaves the 
shell by radiative or conductive means as measured by an observer in the shell and at rest with 
respect to the shells reference frame, which may be written 



Ltot{a + 6a)e 



2{<S>{a+5a)-<S>(a)) 



d<I> 

Ltot{a + 5a) ( 1 + 2— 



+ 2Ltot^ 1 <5a 
da da 



-^^tot(a-) 



^ (T 



(Ltote^*) 



(A3) 



and the energy change during a coordinate time interval, dt, is this times the proper time interval 
e*dt. The two factors of e*(°'+'^")~*(") account for redshift and time dilation, respectively, as the 
energy crosses from the inner to the outer edge of the shell, and we expanded the exponential to 
order 0{5a) as 



„2(<I>(a+5a)-<&(a)) 



d^> 

1 + 2[$(a + 5a) - $(a)l = 1 + 2— 6a. 

da 



(A4) 



Using Eqs. (|Aip . (|A2p and (|A3p in Eq. ([T]) and expressing the change in internal energy in terms 
of the density of internal energy, fjnt, local energy balance can then be written as 



d 



int 



P 



6a 



qe^6adt - Pd ( i(5a ) + adS + d^c 



da 
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(A5) 



Rearranging and noting that 6a is constant in time we then get the gradient of Ltot 
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(A6) 
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Eq. ()A6P is the relation we sought giving the luminosity gradient for stars with variable structure 
and including the contributions from surface and Coulomb terms in the energy density. It may be 
considerably simplified however by noting that we can write the internal energy density in terms 
of the total energy density as S^t /p = e/p — friB in which case from the definition of q 



9 - e t: 



P 



0-* 



de 
dt 



e d/9 
p dt 



(A7) 



dt \ p J p 

Inserting this in Eq. (jA6j) and using the first law as written in Eq. ([3]) finally gives the expression 



P 



de 
di 



Pdp 



p dt 



_a dS 
^ 5a dt 



P 



1 d^c 
6a dt 



a=constant 



ds sr^ dYk 



(A8) 
(A9) 



a=constant 



In this discussion we have used baryon number as the independent variable to emphasize the 
importance of a Lagrangian description for stars with variable structure. This is of course not 
required for all applications and for comparison with other discussions we note that our expressions 
can be converted to use radial coordinate as the independent coordinate through the standard 
relation 

2m' 



da = A-KT p I 1 



dr 



(AlO) 
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